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Climatological  fit  to  the  ionospheric  parameters,  foF2  and  hn,F2  for  the  high 
latitude  stations  at  Sondrestrom,  Greenland  and  Qaanaaq,  Greenland 

1.  Introduction 

Ionospheric  models  such  as  lONCAP  (Lloyd  et  al.  1978),  International  Reference 
Ionosphere  (IRI-90,  Bilitza  1990),  and  PRISM  (Daniell  et  al.  1995),  are  used  for  prediction  of 
ionospheric  parameters.  The  models  lONCAP  and  PRISM  provide  standard  ionospheric 
parameters,  whereas  the  models  IRI-90  and  PRISM  are  used  in  additional  studies  of  the  altitude 
dependence  of  the  ionospheric  electron  density  profile.  Models  such  as  PRISM  provide  a  special 
capability  of  updating  the  ionospheric  parameters  by  using  a  variety  of  real  time  data  collected 
from  various  platforms.  On  the  other  hand,  real  time  operational  systems  need  relatively  simple, 
short  and  robust  algorithms  for  a  prompt  prediction  of  the  ionospheric  parameters  as  a  backup  for 
periods  when  the  system  may  be  temporarily  degraded  or  inoperable.  With  the  availability  of  a 
continuous  chain  of  data  covering  a  complete  solar  cycle,  along  with  high  computer  speed  for 
analyzing  large  volumes  of  data,  it  is  now  possible  to  construct  relatively  simple  algorithms  for 
individual  stations  for  a  prediction  of  ionospheric  parameters  for  the  respective  locations  such  as 
those  needed  for  the  sensor  based  C/NOFS  program  launched  by  the  US  Air  Force.  This  report 
presents  a  construction  of  such  algorithms  for  foF2  and  hniF2  for  the  stations  Sondrestrom, 
Greenland  and  Qaanaaq,  Greenland.  At  the  high  latitude  of  Sondrestrom,  often  a  tongue  of 
ionization  ((TOI)  Tsunoda  1988)  is  seen,  which  is  a  source  for  the  formation  of  polar  cap 
patches.  The  polar  cap  station  at  Qaanaaq  sees  polar  cap  patch  activity  (see  Dandekar  and 
Bullett  1999  and  references  therein,  Dandekar  2002).  Therefore  climatological  algorithms  for 
foF2  to  provide  distributions  corresponding  to  the  presence  or  absence  of  polar  cap  patch  (or  TOI) 
activity  are  also  considered. 

2.  Database 

For  each  of  the  stations:  Sondrestrom,  Greenland  and  Qaanaaq,  Greenland,  four 
representative  one  year  data  sets  of  foF2  and  hmF2  for  the  years  1989- 1990,  1991-1992,  1993-1994 
and  1996-1997  are  used  to  cover  the  whole  range  of  variation  over  solar  cycle  22  (and  the 
beginning  of  solar  cycle  23).  Usually  these  data  were  collected  at  15-minute  intervals  from  the 
digital  ionospheric  sounding  systems  (DISS)  (Reinisch  et  al.  1983)  deployed  at  the  respective 
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stations  by  the  Air  Force  Research  Laboratory,  Hanscom  AFB,  Massachusetts,  US.  The  foF2  and 
hmFi  data  are  extracted  from  the  ionograms  by  the  automatic  scaling  program  ‘ARTIST-4’  of  the 
DISS  system.  (The  accuracy  of  scaling  by  ARTIST-4  will  be  discussed  later).  The  geographic 
and  corrected  geomagnetic  coordinates  of  the  stations,  with  relationships  of  local  time  (LT), 
universal  time  (UT),  and  corrected  geomagnetic  time  (CGMLT)  are  listed  in  Table  1.  Qaanaaq, 
3.3°  away  from  the  corrected  geomagnetic  pole  is  a  real  polar  cap  station,  whereas  Sondrestrom 
at  75.16°  corrected  geomagnetic  latitude  (CGMLAT)  is  an  auroral  station  during  local  day  time 
and  a  polar  cap  station  during  local  nighttime  conditions. 


Table  1.  Stations  used  for  the  study 


Station 

Sondrestrom 

Qaanaaq 

Geogr.  Lat. 

67.0°  N 

77.5°  N 

Geogr.  Long. 

50.6°  W 

68.7°  W 

Corr.  Geomag.  Lat. 

75.16° N 

86.71°  N 

Corr.  Geomag.Long. 

42.47°  E 

40.39°  E 

LT  at  0000  ijT 

2036 

1924 

CGMLT  at  0000  UT 

2154 

2148 

UT  at  1200  LT 

1524 

1636 

Data  availability  from  the  respective  stations  is  presented  in  Table  2.  The  data  selected 
are  typically  from  July  1  to  June  30  of  the  following  year.  In  the  top  two  lines  the  table  lists  the 
station  and  the  parameter  used  for  the  study.  The  left-hand  column  lists  the  period  of  the  data 
used.  For  each  parameter  the  table  presents  the  number  of  observations  for  the  combined  year, 
and  as  percent  available  of  the  expected  data  coverage.  The  table  shows  that  the  data  availability 
was  as  low  as  42%  for  hniF2  at  Sondrestrom  in  89-90  and  as  high  as  90%  for  foF2  at  Qaanaaq  in 
96-97.  The  overall  data  availability  is  62%  for  foF2  and  58%  for  hniF2  at  Sondrestrom,  and  78% 
for  foF2  and  76%  for  hmF2  at  Qaanaaq.  The  averaged  sunspot  numbers  (SSN)  with  the  standard 
error,  for  the  respective  periods  in  Table  2  are  149  ±  49, 127  ±  51,  38  ±  22,  and  10  ±  10 
respectively.  Thus  the  data  used  form  a  good  sample  covering  the  whole  range  of  the  solar  cycle 
change.  Auxiliary  data  such  as  Kp,  sunspot  number  (SSN)  and  solar  flux  (w/cm'^  s'‘,  at  2800 
MHz)  used  in  the  analysis  were  available  from  the  World  Data  Center,  Boulder,  Colorado. 
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Table  2.  Data  available  for  the  study 


Station 

Sondrestrom 

Qaanaaq 

Parameter 

foF2 

hmF2 

foF2 

hmF2 

Year 

No.  of 
obs. 

% 

No.  of 
obs 

% 

No.  of 
obs. 

% 

No.  of 
obs. 

% 

89-90 

15920 

45 

14756 

42 

29698 

85 

29182 

83 

91-92 

21974 

63 

20609 

59 

22545 

64 

21509 

61 

93-94 

26862 

77 

25075 

72 

26308 

74 

25694 

73 

96-97 

22644 

65 

20635 

59 

31548 

90 

30880 

88 

Total 

87400 

62 

81075 

58 

109829 

78 

107265 

76 

3.  Analysis 

3.1  Climatological  fit  to  all  the  data 

As  a  starting  point,  the  determination  of  an  empirical  fit  to  foF2  and  hmF2  data  from 
Sondrestrom  and  Qaanaaq,  without  any  consideration  of  the  TOI  or  polar  cap  patch  activity,  is 
conducted  by  using  an  equation  of  the  form, 
foF2  (or  hniF2 )  =  AO 
+  A1  F 

+  A2sin((H-H0)  *15°*  7t/180“) 

+  A3  cos(0.5  (D-Dl)  (360/366)  *  k  /1 80°) 

+  A4  cos((D-D2)  (360/366)  *  n  /1 80°) 

+  A5Kp  (1) 

where, 

AO  to  A5  are  empirical  coefficients, 

HO,  Dl,  and  D2  are  empirical  constants 
H  is  the  universal  time  (U.  T.)  in  hours, 

D  is  the  incremental  calendar  day,  taking  into  account  the  time  for  the  observation, 

F  is  the  solar  flux  in  w/cm‘^  s'*  measured  at  2800  MHz,  and  averaged  for  90  days  prior  to 
the  day  of  the  ionospheric  observation, 

Kp  is  the  index  of  planetary  magnetic  activity,  and 

the  term  n  /1 80°  converts  the  respective  angles  in  degrees  to  radians. 

The  multiplication  terms  of  coefficients  A1  to  A5  represent  the  dependence  of  foF2  (or 
hmF2)  on  solar  flux,  diurnal,  semiannual  and  annual  variation,  and  dependence  on  the  index  Kp 
of  planetary  magnetic  activity  respectively. 
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All  the  needed  empirical  constants  and  coefficients  in  Eq.(l)  were  determined  by  using 
data  from  all  four  years,  summarized  in  Table  2  with  the  SVDC  routine  of  the  Interactive  Data 
Language  (IDL  1995,  reference  guide)  computer  programs.  For  each  station  and  each  parameter, 
a  contribution  from  every  term  in  Eq.(l)  is  estimated. 

The  terms  contributing  less  than  10  percent  are  dropped,  and  the  remaining  constants  and 
coefficients  are  recomputed.  The  final  terms  used  for  the  fit  are  presented  in  Table  3.  After  all 
the  constants  in  Eq.(l)  are  determined,  one  can  calculate  the  contribution  of  each  term  of  Eq.(l), 
and  its  percent  contribution  to  the  ionospheric  parameter  under  study.  These  relative 
contributions  are  also  presented  in  Table  3.  The  top  two  lines  in  Table  3  list  the  parameter  and 
the  station.  The  next  three  lines  present  the  constants  determined  for  Eq.(l)  for  the  respective 
parameter  and  station.  The  next  six  lines  present  the  computed  values  of  the  coefficients  AO  to 
A5.  For  the  coefficients,  each  column  is  subdivided  into  two  columns.  The  first  column  lists  the 
magnitude  of  the  coefficient  and  the  second  column  gives  the  percent  contribution  of  the 
respective  term  to  that  ionospheric  parameter.  It  is  observed  that  in  Eq.(l),  solar  flux  averaged 
for  90  days  (3  months)  prior  to  the  day  of  ionospheric  observation  provides  a  better  fit  than  that 
based  on  the  use  of  daily  solar  flux.  Therefore  the  coefficients  derived  with  use  of  3  months’ 
averaged  solar  flux  are  listed  in  the  Table  3. 

At  the  bottom  of  Table  3,  each  top  column  is  subdivided  into  three  columns.  The  left- 
hand  column  describes  the  item  for  the  fit  of  Eq.(l).  For  each  parameter  there  are  three  entries; 

1)  standard  error  in  MHz  (or  km),  2)  shift  in  the  fit  with  respect  to  the  observed  values,  and  3) 
standard  error  as  percent.  Similarly  the  errors  between  observations  and  predictions  from  the 
lONCAP  and  PRISM  models  were  computed  and  summarized  in  the  bottom  lines  in  Table  3. 

The  table  shows  that  the  fit  by  use  of  Eq.(l)  (third  line  from  the  bottom  in  Table  3)  has  smaller 
error  than  from  either  the  lONCAP  (second  line  from  the  bottom  in  Table  3)  or  the  PRISM  (last 
line  in  Table  3)  models. 

Correlation  coefficients  were  computed  between  foFa  (or  hmF2)  and  the  respective 
parameter  on  the  right  side  in  Eq.(l)  and  are  presented  in  Table  4.  The  parameters  are  listed  in 
the  left-hand  column.  Two  columns  on  the  right-hand  side  are  for  two  stations.  Each  station 
column  is  divided  into  two  columns  for  the  parameter  foF2  and  hniF2.  The  partial  correlation 
coefficients  are  shown  in  the  table.  Two  lines  at  the  bottom  of  the  table  list  the  total  correlation 
coefficient  and  the  percent  volume  of  the  correlated  data  (C^  to  C,  where  C  is  the  coefficient  of 
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correlation).  The  total  correlation  coefficients  are  very  significant  (>  0.707),  showing  that  the 
correlated  population  ranges  from  43  to  75  percent. 


Table  3.  Empirical  Coefficients  for  Equation  1. 


Parameter 

foF2 

hmF2 

Station 

Sondrestrom 

Qaanaaq 

Sondrestrom 

Qaanaaq 

Constant 

HO 

10.5 

13.25 

8.0 

10.5 

D1 

177 

176 

178 

173 

D2  ^ 

— 

— 

181 

173 

Coefficient 

AO 

2.152 

% 

1.458 

% 

182.39 

% 

199.34 

% 

A1 

36 

0.01658 

40 

0.4798 

24 

0.352 

17 

A2 

1.1056 

32 

0.532 

17 

-29.95 

20 

-20.27 

13 

A3 

14 

1.594 

27 

68.2 

24 

80.6 

28 

A4 

-50.12 

33 

-69.04 

43 

A5 

18  j 

-0.01272 

17 

Standard  error  o 

from 

MHz 

MHz 

% 

MHz 

MHz 

% 

Km 

Km 

% 

Km 

Km 

% 

Obs. 

1.56 

Shift 

mm 

59.6 

Shift 

52.4 

Shift 

Eq.(l)Fit 

1.03 

0.003 

mm 

24 

41.8 

0.04 

18 

37.6 

0.09 

13 

lONCAP 

1.1 

0.16 

26 

1.12 

0.27 

30 

52.7 

35.2 

22 

51.3 

38.2 

20 

PRISM 

1.09 

0.32 

25 

1.03 

0.18 

28 

58.5 

18.1 

23 

46.2 

19.5 

18 

Table  4.  Summary  of  Correlation  Study 


Station 

Sondrestrom 

Qaanaaq 

Parameter 

F0F2 

hn.F2 

foF2 

iim¥2 

Dependence  on 

Partial  coefficient  of  correlation 

Flux 

0.5143 

0.4943 

0.6080 

0.4352 

Diurnal 

0.4799 

-0.3596 

0.2336 

-0.2685 

Semiannual 

0.0791 

-0.2763 

0.2656 

-0.4588 

Annual 

-0.2999 

-0.4925 

Kp 

0.0303 

-0.0153 

Total  Corr.  Coeff. 

0.7504 

0.6599 

0.7441 

0.6972 

Correlated  population  (%) 

56  -75 

43-66 

55-74 

49-70 

Residual  errors  were  examined  to  see  if  these  show  any  systematic  behavior,  which  might 
indicate  a  systematic  time  dependent  under-  or  over-estimate.  Because  no  such  behavior  was 
found,  no  further  correction  to  the  algorithm  could  be  made. 
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3.2  Distribution  of  foF2  with  (or  without)  the  presence  of  polar  cap  patch  (or  TOI)  activity 
As  the  high  latitude  stations  see  the  polar  cap  patch  activity  at  Qaanaaq  and  the  formation 
of  a  tongue  of  ionization  (TOI)  at  Sondrestrom,  one  would  like  to  quantify  foF2  distributions  in 
these  two  groups.  For  the  polar  cap  patch  study,  Dandekar  (2002)  suggested  a  procedure  of 
subtracting  an  averaged  background  for  determining  the  polar  cap  patch  (or  TOI)  activity. 
Dandekar  (2002)  treated  short  duration  peaks  (less  than  30  minutes,  that  is,  at  least  two 
consecutive  observations)  as  jitter  in  the  data.  As  the  high  latitude  station  Qaanaaq  sees  a  lot  of 
short  duration  polar  cap  patches  (seen  from  1 -minute  interval  observations);  the  process  of 
counting  these  peaks  and  their  durations  is  modified  for  this  analysis. 

*  The  goal  here  is  to  separate  the  data  into  two  groups  a)  with  the  presence  and  b)  absence 
of  the  polar  cap  patch  (or  TOI)  activity.  The  monthly  distributions  of  all  foF2  data  are  divided 
into  bins,  each  with  a  width  of  0.1  MHz.  It  is  observed  that  these  distributions  do  not  strictly 
follow  the  Gaussian  distribution  (the  Gaussian  distribution  is  defined  by  two  terms:  the  root 
mean  squared  value  (RMS)  and  the  standard  error  for  the  given  set).  Therefore,  the  mode  and 
mediem  distributions  are  also  used  for  the  analysis.  Because  the  selected  interval  of  0.1  MHz  is 
reasonably  narrow,  and  the  distribution  is  not  necessarily  smooth,  the  modal  (highest  occurrence) 
value  was  determined  by  averaging  foF2  for  the  three  highest  occurrences.  The  values  >  1.1*RMS 
(/mode/median)  foF2  are  assumed  to  be  due  to  the  polar  cap  patch  (or  TOI)  activity.  The  base  of 
the  polar  cap  patch  is  assumed  to  be  equal  to  0.8*  RMS  (/mode/median)  foF2.  Also,  the  duration 
of  a  polar  cap  patch  was  assumed  to  be  <  2  hours.  With  these  criteria,  the  monthly  data  are 
scanned  to  determine  the  distribution  of  foF2  for  the  occurrence  of  polar  cap  patch  activity.  The 
monthly  distribution  of  foF2  in  the  absence  of  polar  cap  patch  (or  TOI)  activity  was  computed  by 
subtracting  this  foF2  distribution  from  the  monthly  distribution  of  all  the  data.  For  the  Gaussian 
distribution,  the  error  of  ±  la  covers  68%  of  the  data  population.  In  using  the  mode  and  median 
methods,  the  corresponding  range  points  were  determined.  As  the  results  from  median  and  mode 
methods  did  not  show  any  systematic  and  significant  differences  from  those  from  the  Gaussian 
distribution,  the  latter  method  (Gaussian  distribution)  is  used  in  this  analysis.  As  an  example  the 
foF2  distributions  for  January  1990  for  Qaanaaq  are  shown  in  Figure  1 .  Three  curves  show  the 
distributions  of  foF2  for  1)  all  data  (dotted),  2)  presence  of  polar  cap  patches  (dashes)  and  3) 
absence  of  the  polar  cap  patches  (dash  dot  dash).  The  corresponding  vertical  lines  mark  the  peaks 
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QAANAAQ,  JANUARY  90 


FOF2  (MHZ) 

Figure  1 .  Distribution  of  foF2  for:  1)  all  data,  2)  with  polar  cap  patch  activity,  3)  in  the  absence  of  polar  cap  patch  activity  for  January 
1990  at  Qaanaaq.  Below  3.8  MHz,  the  curve  for  no  patches  and  the  curve  for  all  data  merge. 


and  widths  for  the  respective  category.  The  curve  corresponding  to  the  presence  of  polar  cap 
patches  is  the  right  hand  part  of  the  curve  and  the  curve  corresponding  to  the  absence  of  the  polar 
cap  patches  is  the  left  hand  part  of  the  curve.  Because  of  the  narrow  width  of  0. 1  MHz  for  the 
bins,  all  these  curves  show  some  jitter.  The  foF2  values  are  also  listed  in  the  figure  for 
convenience.  Three  columns  present  average  -a,  average,  and  average  +a,  respectively,  for  the 
condition  listed  on  the  left  hand  side  of  the  columns.  One  would  expect  higher  foF2  for  polar  cap 
patch  activity  and  lower  foF2  in  the  absence  of  polar  patch  activity  with  respect  to  distribution  of 
all  the  data  as  seen  in  the  figure. 

The  values  of  foF2  for  peak  (RMS)  occurrence  are  determined  for  1)  all  the  data,  2) 
background  (absence  of  activity),  and  3)  the  presence  of  activity.  The  widths  (ct)  for  the 
corresponding  groups  are  determined.  Thus  there  are  two  parameters;  (foF2  peak,  and  width  in 
MHz)  for  each  group.  An  algorithm 

Peak  (or  width  foF2)  =  BO  +  B 1  cos((mon-mon0)  *  FAC  *  n  /ISO”)  +  B2  F  (2) 

is  used, 

where  BO,  Bl,  B2,  are  empirically  determined  coefficients,  and 
monO  and  FAC  are  empirically  determined  constants, 
mon  is  the  month  of  the  observation, 

F  is  the  solar  flux  (averaged  over  90  prior  days)  in  w/cm'^  s  ‘  measured  at  2800 
MHz,  and 

the  term  7t  /1 80°  converts  the  respective  angles  in  degrees  to  radians. 

THE  SVDC  routine  mentioned  above  is  used  to  determine  the  empirical  constants  in 
Eq.(2).  For  simplicity  monO  =  6,  and  FAC=15  is  used  for  the  whole  data  set.  The  results  are 
summarized  in  Table  5. 

4.  Interpretation: 

The  distance  between  the  stations  Sondrestrom  and  Qaanaaq  is  1306  km.  The 
asymmetrical  characteristics  of  the  auroral  oval  with  respect  to  the  magnetic  north  pole  makes 
Sondrestrom,  at  75.16°  corrected  geomagnetic  latitude  (CGMLAT),  an  auroral  station  during 
local  daytime  and  a  polar  cap  station  during  the  local  nighttime  conditions,  whereas  Qaanaaq,  at 
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Table  5.  Empirical  Coefficients  for  Equation  2 


Station 

Sondrestrom 

Qaanaaq 

Coefficient 

Item 

BO 

B1 

B2 

Standard 

Error 

(MHz) 

B1 

B2 

Standard 

Error 

(MHz) 

Peak 

F0F2 

(MHz) 

All 

Data 

2.088 

0.855 

0.015 

0.308 

1.301 

1.590 

0.016 

0.423 

Patches 

2.789 

0.344 

0.019 

0.426 

1.731 

1.243 

0.019 

0.486 

Back- 

Ground 

1.766 

1.14 

0.013 

0.265 

0.013 

0.388 

Width 

(MHz) 

1 

All 

Data 

0.628 

-0.364 

0.006 

-0.459 

0.005 

0.198 

Patches 

0.616 

-0.515 

-0.366 

0.004 

0.164 

Back- 

Ground 

0.374 

-0.22 

0.004 

0.214 

3.3°  away  from  the  corrected  geomagnetic  pole,  is  a  true  polar  station  at  all  times.  These  two 
stations  see  different  ionospheric  features.  Sondrestrom  sees  the  tongue  of  ionization  (TOI). 
Under  favorable  interplanetary  magnetic  field  (IMF)  conditions,  the  TOI  breaks  into  polar  cap 
patches  and  these  patches  drift  northward  (poleward)  to  Qaanaaq  (see  Dandekar  and  Bullett 
1999,  Basu  and  Valladares  1999,  and  references  therein).  Against  this  background  let  us  look  at 
the  results  summarized  in  Table  3. 

The  peak  in  diurnal  variation  of  foF2  (or  hn,F2)  is  a  quarter  of  the  day;  that  is  six  hours 
after  HO,  and  is  to  be  added  to  the  constant  HO.  Thus  according  to  Eq.(l)  foF2  peaks  at 
Sondrestrom  around  1630  (10.5+6)  UT,  and  the  peak  at  Qaanaaq  occurs  around  1915  (13.25+6) 
UT.  The  delay  in  the  occurrence  of  the  peak  of  foF2  between  the  two  stations  is  consistent  with 
the  observation  that  the  polar  cap  patches  are  formed  from  the  tongue  of  ionization  (TOI)  at  the 
latitude  of  Sondrestrom,  and  drift  towards  the  pole  to  Qaanaaq. 

The  constants  D1  and  D2  for  the  semiannual  and  annual  variation  are  practically  the 
same.  Note  that  the  annual  component  is  weak  (absent,  A4~0)  for  foF2  at  both  the  stations.  The 
constant  A1  is  same  for  both  the  stations.  The  presence  of  the  constant  A5  shows  a  magnetic 
control  (Kp)  on  foF2  at  both  stations. 

The  respective  coefficients  for  hmF2  for  both  the  stations  are  in  reasonable  agreement 
with  each  other.  Note  that  altitude  dependence  shows  an  annual  component  (A4),  but  weak 
magnetic  dependence  (A5~0). 
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The  bottom  section  of  the  table  compares  the  fit  of  Eq.(l)  with  that  from  the  lONCAP 
and  PRISM  models.  The  lONCAP  and  PRISM  models  are  based  on  data  mostly  collected 
during  solar  cycle  19  (International  Geophysical  year  IGY  1957).  The  table  shows  that  the 
standard  error  (both  in  magnitude  and  as  percent  error)  as  well  as  the  shift  is  smallest  from  the 
algorithm.  The  lONCAP  model  on  the  average  underestimates  the  altitude  by  35  km.  The 
PRISM  model  underestimates  it  by  15  km,  in  spite  of  the  fact  that  it  incorporates  quite  a  few 
more  parameters  than  those  used  in  the  lONCAP  program  for  ionospheric  predictions.  The 
strength  of  the  PRISM  model  lies  in  the  fact  that  the  results  can  be  updated  by  using  real  time 
data  from  various  platforms. 

A  correlation  study  was  conducted  between  the  items  (multiplying  terms  for  the 
coefficients)  listed  in  Table  3  and  the  observations  of  foF2  and  hniF2.  The  results  are  summarized 
in  Table  4.  The  table  presents  partial  correlation  coefficients  for  the  items  shown  in  the  left-hand 
column  for  the  stations  listed  on  the  top  line,  and  for  the  parameter  listed  on  the  second  line  from 
the  top.  The  bottom  line  shows  the  total  correlation  coefficient  for  the  parameter  listed  on  line  2 
from  the  top,  with  the  relevant  items  in  the  left-hand  column  used.  Note  that  items  contributing 
less  than  10%  (see  Table  3)  are  not  used.  The  larger  the  magnitude  of  a  partial  correlation 
coefficient  in  a  given  column,  the  larger  the  contribution  the  item  makes  to  the  parameter.  The 
second  row  from  the  bottom  lists  the  correlation  between  observations  and  the  fit  from  Eq.(l). 
This  row  shows  that  when  the  most  relevant  items  are  combined,  they  produce  a  significant  level 
of  correlation  (>0.707),  indicating  an  overall  improvement  in  the  dependence  of  the  parameter  on 
the  chosen  terms.  The  last  row  presents  the  range  (C^  to  C,  where  C  is  the  coefficient  of 
correlation)  of  percent  population  showing  covariation.  At  Sondrestrom  this  range  is  a  little 
lower  for  hmF2.  All  others  have  a  good  range  (at  least  50%)  of  the  population  showing  a  good 
co-variation. 

Equation  2  is  used  to  determine  the  coefficients  presenting  the  distribution  of  general 
population  and  background  (absence  of  polar  cap  patches)  and  polar  cap  patch  activity  at 
Qaanaaq  and  TOI  activity  at  Sondrestrom.  These  are  summarized  in  Table  5.  The  table  is 
divided  into  three  major  columns,  the  first  for  listing  the  items  and  the  other  two  for  the  stations 
Sondrestrom  and  Qaanaaq.  The  first  column  is  subdivided  in  two  columns;  the  first  lists  the 
parameter  under  study  and  the  second  column  describes  the  category.  For  each  station  the  major 
column  is  subdivided  into  four  columns  to  present  the  empirical  coefficients  BO,  Bl,  B2,  and  the 
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standard  error  for  the  empirical  fit  to  the  data.  For  simplicity  the  reference  month,  monO,  and  the 
factor  FAC  in  Eq.(2)  are  fixed  at  6  and  15  respectively  for  all  the  data  in  the  set.  The  term 
multiplying  the  coefficient  B2  is  F,  the  solar  flux.  As  the  90-day  average  solar  flux  before  the 
day  of  ionospheric  observation  is  found  to  produce  better  results  than  that  from  using  daily  solar 
flux,  the  former  is  used  in  the  computations. 

The  results  are  shown  in  Figures  2A  and  2B  for  Sondrestrom  and  Qaanaaq  respectively. 
Each  figure  has  six  panels.  Each  row  presents  results  for  all  data,  for  polar  cap  patch  activity  and 
absence  of  polar  cap  patch  activity.  Along  the  column  are  the  parameters:  peak,  and  width  for 
the  category  of  the  row.  In  each  panel  the  horizontal  axis  presents  months  of  observation  as  a 
continuous  sequence,  and  the  numbers  just  above  the  axis  show  the  calendar  year  for  the 
observations.  In  each  panel  the  continuous  curve  shows  the  observed  dependence  and  the  dotted 
curve  shows  the  empirical  fit.  In  each  panel  the  numbers  on  the  left-hand  side,  from  the 
top  to  the  bottom,  list  the  coefficients  and  standard  error  respectively  in  the  empirical  fit  of 
Eq.(2).  In  each  panel,  the  comparison  of  the  pairs  of  curves  shows  that  the  empirical  fit  is  very 
reasonable,  in  spite  of  a  factor  of  2  variation  in  both  the  peak  and  the  width  of  the  distribution  of 
foF2  for  the  respective  category,  over  a  period  of  9  years.  Thus  these  algorithms  would  be  very 
useful  to  operational  systems  for  estimating  foFa  when  real  time  observational  data  are  not 
available. 

In  the  above  analysis,  the  automatically  scaled  data  from  the  DISS  systems  are  used 
under  the  assumption  that  this  scaling  is  reasonably  accurate.  For  validating  this  hypothesis,  the 
auto-scaled  data  were  compared  with  manually  scaled  data.  The  results  are  presented  in  Figure 
3.  The  four  panels  in  Figure  3  present  observed  errors  as  histograms  for  foFa  and  hmF2  for  the 
stations  Qaanaaq  and  Sondrestrom.  The  data  for  one  week  of  November  and  one  week  of 
December  1993  were  used.  The  selected  year  is  listed  at  the  top  of  each  block.  The  left-hand 
top  comer  in  each  block  lists  the  number  of  ionograms  available  for  comparison.  The  results  are 
summarized  in  Table  6.  The  table  shows  that  accuracy  in  auto  scaling  of  frequency  is 
acceptable,  whereas  the  accuracy  in  auto  scaling  of  altitude  is  marginal.  These  errors  would 
superpose  on  the  errors  in  the  algorithm  shown  in  Tables  3  and  5. 
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Table  6.  Population  of  data  scaled  using  the  ARTIST-4  program  for  a  given  range  of  error 


Parameter 

Error  (MHz) 

foF2 

Error 

(km) 

hmF2 

Station 

Sondrestrom 

Qaanaaq 

Sondrestrom 

Qaanaaq 

Year 

1993 

1989 

1993 

1989 

Percent  data  covered 

Percent  data  covered 

0.5 

81 

85 

20 

30 

33 

1.0 

93 

91 

40 

48 

52 

1.5 

96 

95 

60 

61 

66 

80 

77 

100 

77 

81 

5.  Conclusion 

Customized  algorithms  derived  by  empirical  fit  to  historical  data  for  individual  stations 
work  better  than  the  models  used  for  global  predictions.  The  next  possible  step  in  improving  the 
present  algorithm  in  Eq.(l)  is  to  subdivide  the  diurnal,  semiannual,  and  annual  variations 
presented  by  a  single  (sine,  cosine)  function,  into  several  sections  by  using  a  spline  technique 
(Boor  1978).  Another  possibility  is  to  put  additional  parameters  into  the  empirical  equations  to 
improve  the  accuracy  of  prediction.  These  relatively  simple  algorithms  should  be  incorporated 
in  models  like  PRISM  for  improving  global  predictions. 

The  high  latitude  stations  selected  here  also  see  special  features  such  as  the  tongue  of 
ionization  at  Sondrestrom  and  polar  cap  patch  activity  at  Qaanaaq.  An  additional  algorithm  in 
Eq.(2)  quantifies  the  foF2  distributions  for  the  respective  periods,  and  would  be  useful  to 
operational  systems,  which  are  often  degraded  during  disturbances  associated  with  these 
phenomena. 
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